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Abstract. This article settles the convergence question for multivariate bary- 
centric subdivision schemes with nonnegative masks on complete metric spaces 
of nonpositive Alexandrov curvature, also known as Hadamard spaces. We es- 
tablish a link between these types of refinement algorithms and the theory of 
Markov chains by characterizing barycentric subdivision schemes as nonlinear 
Markov semigroups. Exploiting this connection, we subsequently prove that 
any such scheme converges on arbitrary Hadamard spaces if and only if it con- 
verges for real valued input data. Moreover, we generalize a characterization 
of convergence from the linear theory, and consider approximation qualities 
of barycentric subdivision schemes. A concluding section addresses the re- 
lationship between the convergence properties of a scheme and its so-called 
characteristic Markov chain. 



Introduction and main results 

The convergence and smoothness analysis of refinement schemes processing data 
from manifolds and more generally metric spaces has been a subject of intense 
research over the last few years. As to convergence, complete spaces of nonpositive 
curvature have proven most accessible in terms of generalizing well-known facts 
from the linear theory to the nonlinear setting. An example of such a structure 
prominent in applications is the space of positive definite symmetric matrices, which 
represent measurements in diffusion tensor imaging. 

While the question whether the smoothness properties of the linear model scheme 
prevail when passing to the nonlinear setting was successfully addressed in [6], the 
corresponding convergence problem remained unsolved. The aim of this article is 
to fill this gap in the theory. 

Relying on a martingale theory for discrete-time stochastic processes with val- 
ues in negatively curved spaces developed in [TU], we observe that the refinement 
processes in question actually act on bounded input data as nonlinear Markov 
semigroups. This fact will substantially facilitate their convergence analysis. 

Let us specify the general setup. Given a metric space {X,d), a refinement 
scheme is a map S : ^°°(Z",X) -> ^°°(Z^X). We caU S convergent if for all 
X G £°^{Z'',X) there exists a continuous function S°"x -.W^X such that 

doo{S°^x{-/2"),S"x) =snpd{S'^x{j/2"),S"xj) -^0 as n ^ oo. (1) 

j 

Visualizing S" function on the refined grid 2 "Z^, convergence to S°°x is 

tantamount to <ioo('5'"a;, /|2-iz0 ~^ 0. 
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Throughout the present paper we are mostly concerned with so-caUed barycen- 
tric refinement schemes associated to nonnegative real-valued s-variate sequences 
(ai)iez= of finite support, henceforth referred to as masks, which we require to fulfill 
the basic sum rule 

a,_2i = 1 for i G Z'*. (2) 

Barycentric refinement schemes act on data x G £°°{'Z^,X) from a complete metric 
space of nonpositive curvature in the sense of A. D. Alcxandrov according to the 
following rule: 

Sxi = argmin ^ ai^2j(f{xj, •) • (3) 
Vez= J 
Much is known about the convergence of these type of refinement algorithms in 
the case AT = R. On complete manifolds of nonpositive sectional curvature conver- 
gence analysis was initiated in the article |12| . The author in [5j recently proved 
general convergence statements for arbitrary Hadamard spaces using the principle 
of contractivity: A scheme is called contractive with respect to some nonnegative 
function I? : £°°(Z*, A") IR+ if and only if there is 7 < 1 such that 

D{Sx) <-fD{x), for all x££°°{Z\X). 

The function D is referred to as a contractivity function for S. An important class 
of contractivity functions is associated to balanced, convex and bounded subsets fl 
of R*: 



Dn{x) = sup d{xi,Xj), (4) 

p(i-j)<2 

where p denotes the Minkowski functional of fl. Contractivity functions of this type 
are called admissible, cf. [5]. The following result is taken from loc. cit.: 

Proposition 1. A barycentric refinement scheme with nonnegative mask which is 
contractive with respect to some admissible contractivity function also converges. 
This implies convergence in case the support of the mask coincides with the set of 
lattice points within a centered unimodular zonotope or a lattice quad with nonempty 
interior. 

The main result of the present article is a substantial extension of this statement 
and describes a phenomenon which could be referred to as linear equivalence: 

Theorem 1. A barycentric refinement scheme converges on arbitrary Hadamard 
spaces if and only if it converges on the real line. 

The proof of this fact, given in Section [51 relies on a stochastic interpretation of 
the subdivision rule ([3]). More precisely, for each nonnegative mask a = (ai)iezs 
satisfying the basic sum rule ^ one finds a so-called characteristic Markov chain 
with state space 7/ and transition matrix (ai^2j) jjez= in terms of which the 
iterates of the refinement algorithm acting on x G ^°°(Z*, A") may be written as 

S''x,^E(xoXl\\\X^^i), 

see Theorem m Here E{ - |||A^o) denotes the filtered conditional expectation intro- 
duced by K.-T. Sturm in [10]. Thus, as in the linear case. 

No ->Lipl(£-(Z^A)); 



NONLINEAR MARKOV SEMIGROUPS AND REFINEMENT SCHEMES 



3 



may be considered a (nonlinear) Markov semigroup. Here Lipi{£°° {Z^ , X)) refers 
to the set of maps T : £°°{Z'^,X) — > £°"{'Z'^,X) satisfying the Lipschitz condition 

d^{Tx,Ty) <d^ix,y) for x, y £ £°°(Z^ X). (5) 

Combining Theorem [T] with other recent developments in the theory of linear sub- 
division schemes with nonnegative masks and their barycentric counterparts on 
nonlinear objects, one comes up with a variety of remarkable results: 

In the articles [ini and [13] , X. Zhou establishes general theorems on the relation 
of the mask's support with its convergence properties, which, utilizing Theorem [T] 
now generalize to the following: 

Theorem 2. A barycentric subdivision scheme S : £°°{Z^,X) — !> £°°{'E'^,X) with 
nonnegative mask converges under each of the following circumstances: 

(i) The support of (ai)i^i' coincides with the set of grid points inside a balanced 
zonotope. 

(ii) The grid dimension s = 1 and, if, after a possible index translation, {ai)i^z = 
(. . . , 0, 0, Oq, . . . , Cat, 0, 0, . . . ), the integers within the support are relatively prime 
and < ao,ai < 1. This also constitutes a necessary condition for convergence. 

Moreover, as far as finite-dimensional Hadamard manifolds are concerned, the 
smoothness question is settled by a combination of Theorem [T] and recent work 
from [6]: 

Theorem 3. On a smooth Hadamard manifold, a barycentric subdivision scheme 
S : £°°{'Z,'^,X) £°°{lf,X) with nonnegative mask converges and produces r -times 
differentiable limit functions if and only if the same is true for the corresponding 
linear scheme. 

1. Stochastic preliminaries 

This section is devoted to a stochastic interpretation of the subdivision rule 
More precisely, we view barycentric subdivision as the semigroup acting on 
£°°('Zi'',X) associated to the so-called characteristic Markov chain of (ai)igz=- This 
result requires some prerequisites about the notion of conditional expectation of 
random variables with values in Hadamard spaces. This type of metric space was 
first investigated by A.D. Alexandrov in his seminal articles [IJ and [2]. In loc. cit., 
Alexandrov uses the Gauss-Bonnet theorem, i.e. the fact that the deficiency of the 
angle sum in a geodesic triangle can be expressed by means of the ambient space's 
curvature, to generalize the notion of curvature bounds to a purely metric setting. A 
comprehensive introduction to the nowadays well-established theory of such spaces 
is [3]. Furthermore we refer to the survey article [llj on probability measures and 
their centers of mass on Hadamard spaces. A treatise of the smooth case can 
be found in [7i, which investigates infinite-dimensional Riemannian manifolds of 
nonpositive curvature. 

Definition 1 (Hadamard spaces). A complete metric space is called Hadamard- or 
global NPC-space if and only if for xo,xi one finds y ^ X such that the so-called 
Hadamard inequality holds true for all z G X: 

d{z, yf < ^d{z, xof + ^d{z, xif - ^d{xn,xif. (6) 

In other words, a Hadamard space is a complete metric space with nonpositive 
curvature in the sense of A. D. Alexandrov. The Hadamard inequality ^ expresses 
the fact that geodesic triangles are 'slim' compared to Euclidean ones with the same 
edge lengths. 
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Recall that a continuous curve xt & X , < t < l,is called geodesic if and only if 
d{xs, Xt) — d(xo, xi)\t ~ s\ for all < s, t < 1- A well-known property of Hadamard 
spaces is that they are strongly geodesic: given any two points xo,xi S X, one 
finds a unique geodesic joining them. In particular, there is meaningful notion 
of convexity. More precisely, a subset K C X oi a global NPC space is called 
convex if and only if for each two xo,xi G K, the joining geodesic xt remains in K. 
Nonpositive curvature turns out to have a major impact on the convexity properties 
of the metric d( ■ , ■ ). As announced above, we are particularly interested in the 
concept of conditional expectation for random variables with values in Hadamard 
spaces. Although several approaches appear in the literature, a definition due to 
K.-T. Sturm, see [lOj, serves our purposes best. Sturm's definition relies on a convex 
projection property enjoyed by NPC spaces: 

Proposition 2. Suppose {X, d) is a global NPC space, and C C X a closed and 
convex subset. Then there exists a well-defined projection map ttc : X — > C defined 
by the relation 

d{Trc{x), x) — mm d{y , x) . 

This projection is Lipschitz- continuous in a sense that d{TTc{x),'Kc{y)) < d(x,y) 
forx,y € X. 

Given a probability space (ri,ij, P) and a metric space {X,d), a strongly mea- 
surable function is called square-integrable if and only if 

[ d^{f{uj),x)P{duj) < oo 

for one (and then all) x E X. The space L'^{^,X) of square-integrable functions 
f : ^ X factorized by the relation of being equal almost surely is endowed with 
a metric 

In case {X,d) is a global NPC space, it is well-known that X), ^2) inherits 

the Hadamard property. Moreover, given a subalgebra 5 C 3^, it is easy to show 
that L^{Q,X), the subspace of L'^{^,X) corresponding to ^/-measurable square- 
integrable maps, is closed and convex. In view of Proposition [2] one thus obtains 
the following Definition: 

Definition 2 ([10]). Suppose is a probability space, and let Y e L^{^,X) 

be a square-integrable random variable with separable image in the Hadamard space 
X. Given a subalgebra Q Q ^, there is a ^/-measurable Z . ^ X minimizing the 
functional 

Z' ^ E{d^{Y,Z')) 

among all ^/-measurable square-integrable random variables. Any other minimizer 
coincides with Z almost surely. One refers to Z as the conditional expectation of Y 
given Q and uses the notation Z = E{Y\Q). 

Otherwise put, following the notation of Proposition [21 E{Y\Q) = 7ri2(g x)(^)- 
This definition follows the principle that in the real-valued case, the conditional 
expectation E{-\g) : Li(5^,]R) L'^{g,R) as introduced by Kolmo gorov restricts 
to as the orthogonal projection to the space of C/-measurable L^-functions. 
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Remark 3. The space LP{^,X) is, as in the case p — 2, defined to be the space 
of strongly measurable functions / : (r2,5^,P) X such that J^dP{f{uj),x)F{duj) 
modulo equality almost surely. Again LP{^,X) comes with a metric of the form 

The conditional expectation as defined above is continuous in a sense that for each 
two Y,Z e L'^i'S, X) and p e [1, oo], 

dpiEiY\g),E{Z\g)) < dp{Y,Z). 

In particular, E{ ■ \Q) extends continuously to L^{'^,X). 

Remark 4. It is inherent in the existence statement provided by Definition [5] that 
each integrable random variable with values in X possesses an expected value defined 

by 

E{Y) :=i?(y |{0,r!}), 
which is the unique minimizcr of the functional 

( d^{Y{Lu),z)F{dijj). 
Jn 

Recall the elementary smoothing property of conditional expectations for random 
variables with values in a linear space: Given an ordered pair of subalgebras Q C 
(y C 3^, one has 

EiEi-\g)\g) = Ei-\g) 

almost surely. Not surprisingly, this associativity property is violated in the non- 
linear case, see Example 3.2 in [TOj. Actually a generalization of this feature would 
render the theory of nonlinear subdivision schemes of type ^ obsolete, as the 
discussion following Theorem [S] below illustrates. 

It is for the reason of lacking associativity that K.-T. Sturm in his treatise |10| 
defines: 

Definition 3. Suppose 5o ^ S^i ^ • ■ -Sn = 5^ is a sequence of subalgebras. Fur- 
thermore assume Y € L'^{^,X). Then one defines the filtered conditional expecta- 
tion of Y given {^n)o<n<N as 

E{Y\\\do) =E{--- E{E{Y\^N-i)\dN-2) ■ ■ ■ |5o). (7) 

Let us briefly recall the basic construction of a Markov chain, beginning with a 
foundational definition. Recall that a family of maps 

Pn^rn : X ^ R 

parametrized by nonnegative integers n < m, is called Markov transition kernel if 
the following are fulfilled: 

(i) Pn.m{i, ■) is a probability measure for all i 
(u) p„^„(i, •) = (5{i}. 

(iii) The Chapman- Kolmogorov equations are fulfilled: for each n < £ < m one 
has 

) ^ ^ Pn,S^k)pLrn{k,j)- (8) 

We use the notation P — {pn.m)n<m- 

Let V, = (Z*)'*'° denote the space of sequences on the grid Z'*, endowed with the 
infinite power ^ = 0^gpj*P(Z'*) of the discrete sigma algebra on Z". Consider the 
filtration 

= ct({zo} X • • • X {i„} X Z'' X • • • I e Z" for j = 1, . . . , n). 
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Choose an initial distribution a awlf , and introduce a probability measure on 
via 

and the standard extension theorems. Whenever a equals the point measure 
on {z}, we will write Pq, = P^. Expected values of integrable random variables 
y : n ^ X with respect to Pq are, as usual, written as 

Finally, the discrete stochastic process 

constitutes a Markov chain associated to the transition kernel P, meaning that the 
linear Markov property holds true; For any nonnegative f : 1/ one has 

i?a(/(^rn)|5n) = ^ Pn,™ (^n , J )/ (j ) • (9) 

In particular, this linear Markov property allows for the interpretation of Pn,m{hj) 
as the transition probability Va{Xm = j\^n — i)- 

In view of the convergence analysis of barycentric subdivision schemes, it is 
of particular interest to gain a deeper understanding of principle of conditioning 
in case the filtration stems from a Markov chain. A nonlinear Markov property 
analogous to ([§]), see [101 Theorem 5.2], leads to a representation of the conditional 
expectation explicit enough for our purposes. We provide a short proof adapted to 
our setting, beginning with an auxiliary result which can be found e.g. in jlUj : 

Lemma 5. Suppose (Xfe)fcgN(j is a Markov chain in 'L^ associated to the transition 
kernel P. Choose an initial distribution a. Furthermore assume Y : ^ X 
is ^n-fneoLSurahle, and let x G i°°{'Z/^X). Then for nonnegative and measurable 
/ : X X X — > M and m > n we have 

[ f{xoXr,,{io),Y{uj))F^{dLo)= [ Vp„,™(X„H,j)/(2;,,rH)P«(dw). 
Jn Jn J 

Proposition 6 (Nonlinear Markov property). Let {Xk)ket>!o ^ Markov chain as 
in Lemma\^ and suppose x G ^°°(Z*, X), with {X, d) Hadamard. Choose n,m gNq 
with n < m. Then 

Ea{x{Xjn)\dn){i^) = argmiu ^ p„,„(X„(w), j)d^(xj, • ) 

jez= (10) 

= Ex^{u){x{Xjn))- 

Proof. By the linear Markov property © , 

Y{lj) := argmin {E„{d^{x o ■ ) | 5„)(a;)) 
= argmin ^ p„,,„(X„(cj), j)d^(xj , • ). 

Clearly Y , as a measurable function of X„, is 5^„-measurable. Thus, in order to 
verify that Y is indeed the conditional expectation of x{Xm) given it remains 
to show that for each Jn-nieasurable function Z the inequality Ea{d^{XrmY)) < 
Ea{d'^{Xra, Z)) holds truc, cf. Definition [21 For this sake, define 

V' : Z" X X ^ R>o U {oo}; 
{i,z) i-^J2jPnMi,j)d{xj,z). 
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By construction of Y we have tp{X„i,Y) < ip{X„i, Z). Thus, Lemma [S] implies 

Ejd\xoX^,Y))=E^{,l;{X„,,Y)) 
<E^i^P{Xm,Z)) 
= E^{<f{xoX,^,Z)). 

□ 

Remark 7. Proposition [5] imphes that the expression Ea{x{Xm)\'Sn) actuaUy is 
independent of the initial distribution a. Therefore we will omit a in the following 
and simply write E{x{Xm)\^n)- 

We are now in a position to establish a link between nonlinear Markov semigroups 
and barycentric refinement processes. Suppose a = (ai)igz= is a nonnegative com- 
pactly supported s-variate sequence such that ai-2j = 1 for all « G Z*. Define 
recursively = (5{o} (*) E^nd 

Then clearly 

a r ■ ■\ (m—n) 

defines a Markov transition kernel. This kernel is homogeneous in the sense that 
Pn,n+i{hj) = ai-2j for any n E Nq. We write P'^ = (Ki,m)n<m, denote the 
associated Markov chain by X^, and refer to X^ as the characteristic Markov chain 
for (ai)igz=- The central observation of this article is the following consequence of 
the nonlinear Markov property PTI| : 

Theorem 4. Suppose x : Z,'^ X is bounded, where {X, d) is a Hadamard space. 
Let S be a barycentric refinement scheme acting on data from X according to the 
subdivision rule Q. Let X^ denote the characteristic Markov chain of {ai)is^is . 
Then 

S^^xoX^^EixoX^mo). 

Proof. This statement is proven by induction over n using the following computa- 
tion and Proposition [BJ 

EiS'^-'^x o X^ I = argmin ( ^ pk-i.k{Xk-i,j)d\S^~''x„ •)) 

= argmin ( ^ axj_,-2jrf'(^""''2;, , •)) 
= S''-''+^xoX^^^. 

□ 

The remainder of this article is devoted to analyzing the effects of this repre- 
sentation of the iterates of S on the convergence properties of barycentric schemes 
with nonnegative masks. 

2. The convergence problem 

We begin this section by summarizing some well-known facts about the conver- 
gence of barycentric schemes acting on real-valued input data. Classical resources 
on this topic are [51 HI [S] . 
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Theorem 5. Suppose a = (ai)igzs is an s-variate compactly supported sequence of 
nonnegative reals. Define a refinement scheme S : £°°(Z*,M) — >■ i?°°(Z*,R) via 

Sxi — ai^2jXj, where x £ ^°°(Z*). 

Then a necessary condition for the convergence of S onM. is the basic sum rule ^ . 
In case the mask {oiji^z' obeys this rule, we conclude 

Sxi = argmin ^ ai-2j\xj ^ ' M = argmin ^ ai^2jdy\{xj, • )^ 

Moreover, S converges if and only if there exists a continuous : M — > M subject 
to the functional equations 



v{t) = Y.^'M'^t-j) (11) 
E^(^-j) = 1- (12) 

i 

Due to Equation pip . ip is referred to as an SL-refinable function. Given bounded, 
real-valued input data {xi)ii£z^, the limit function may be written as 

S°-x{t)=Y,^{t~j)x,. 

In particular, ip = S°°d^Qy, where (5{o} denotes the Dirac distribution on the origin. 

Assuming that conditional expectations of bounded random variables mapping 
to the metric space X are well-defined in the sense of Definition [2j and in addition 
satisfy the smoothing property ([7]), we could deduce from Theorem 3] 

S"xoXo^E{xoX^mo) 
^EixoX^^l^o) 

= argmin(^p(")(Xo,j)d'(a^:;-, •))• 
j 

Note, however, that pf^{i,j), the n-step transition probabilities of (X^)/jgNoj can 
be viewed as (S'"(5o)i_2"j, where S denotes the linear counterpart to S, and 
the Dirac delta on the origin, cf. Theorem [SJ Thus, the assumption of associativity 
would guarantee that every scheme converging for linear input data would converge 
on X as well. Indeed, the limit functions for given input data x £ £°°{Z'' ,X) would 
satisfy 

S^x{t) = aTgmin{J2^{t - j)d'{x„ ■)), 

3 

where Lp = ^""^{o}, leading to a complete analogy to the linear case. However, non- 
linear conditioning is a non-associative operation. The above observations demon- 
strate that this lack of property ([7]) constitutes the need for a further discussion of 
convergence. 

The first result of this section is a small, but useful generalization of Theorem 
1 in [S]. Although the proof transcribes more or less directly, we give a detailed 
exposition for the reader's convenience. 
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Theorem 6. Let S, T be refinement schemes acting on data from a Hadamard 
space X . Then S converges under the following assumptions: 

(i) There is a function D : t°°[l/ ,X) — > R>0; a nonnegative real number 7 < 1 
and a positive integer no such that 

DiS'^x) <-f^"/''°^D{x) (13) 

forx e e°°{Z^,X) and n e N. 

(ii) T is convergent and satisfies 

{Tx,Ty) < doo{x,y) 

forx,y e ^°°(Z^X). 

(Hi) There is C > such that 

doo{Sx,Tx) <C -Dix) 

forxe £°°(Z^J^). 

Proof. We set := T°°(S'"a;)(2"y) and claim that this defines a Cauchy se- 

quence in {C{R^,X),doo)- Note first that given n S N and ?/ e M'', by continuity 
of /„ respectively /„+i, we find j G and m G N such that 

d(/,(2/),/,(2-"j)) <C-I?(a;)7["/"°1 iorr = n,n + l. (14) 

Moreover, due to convergence of T, by multiplying both the numerator and the 
denominator of the number with a power of two if necessary we may assume 

m to be sufficiently large for 

d(/r(2""j),r™-''(S""a;)j) = d(T°°S"^a;(2'-™j),r™-''(S"'x)j) < C ■ D{x)j^''/''°^ 

to hold for r = n, n + 1, in addition to (fH|) . This together with (i) and (ii) implies 

d(/„(y),/„+i(2/)) < d(/„(y),/„(2-'"j)) 

+ d(/„(2-™j),r™-"5"a;,) 

+ d{T"'-"S''xj,T'"-'''-^S"+^Xj) 

+ d(r™-"-i5"+i.x,,/„+i(2-"j)) 

+ d(/„+i(2-™j),/„+i(2;)) 

< 5C-i:>(a;)7["/""l, 

showing that /„ is a Cauchy sequence. Since X is complete, we find a continuous 
f -.W ^ X with fn^f uniformly. We claim that S"x converges to / in the sense 
of (dJ. For m>n and j G Z*, we obtain the inequality 

m—l 

d{T"'-"S''xj,S"'xj) < J2 diT'^'^^S^Xj.T'^'-^-^S^+^j) 

k—n 

which together with 

d(/„(2-™j),5"x,) < d(/„(2-"j),T'"-"5"x,) + d(T"-"^"a;j,5™a;,) 
establishes the claim. □ 

Definition 4. In accordance with [S], we call a scheme S satisfying weakly 
contractive. Thus, a weakly contractive scheme is contractive if and only if tiq = 1. 
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In the following we rely on a nonlinear version of Jensen's inequality due to K. 
T. Sturm, of. [TO]: 

Lemma 8 (Conditional Jensen's inequality). Suppose ^ : X ^ M. is a convex, lower 
semicontinuous function on a Hadamard space {X,d), and (17,5^, P) is a probability 
space. Moreover suppose {^k)keNo filtration in ^. Then for each bounded, 

^ N -measurable random variable Y : Q ~^ X the following holds true: 

^{E{Y\mk)k>n)) < E{i,{Y)\Sn). (15) 

Lemma 9. Suppose the linear scheme associated to (ai)i,=z= converges, a7idsupp(a) 
C fl, with Q bounded, convex and balanced. Denote by p : R'^ ^ M>o the Minkowski 
functional of CI. Furthermore define D : £°°(Z,'^,X) — > M>o via 

D{x) = sup d{xi,Xj). 

p{i-j)<2 

Then the barycentric scheme S associated to (ai)igzs *s weakly contractive with 
respect to D. 

Proof. The Hadamard property implies that for each zq £ X the function 

X — > M>o; zi— >d(z,zo), 

which clearly is continuous, is convex as well. Thus, by Jensen's inequality (llSp 
and Theorem m 

d{S^xoXo,zo)^d{E{xoX^Mo),zo) <Eid{xoX:^,zomo)- (16) 
Recall that the transition kernel of X^ takes the form 

\ / n<m 

Thus Proposition [S] implies i?(d(X^, zo)|S^o) = X^feez^ '^jf«-2"fe'^(^fe' -^o)- Together 
with this gives 

d{S'^Xi,ZQ) < ^ a-"2r.A;d(a;fe, zo) for all i G Z''. 

Substituting S^^Xj for zq, we deduce 



d{S"x,,S''xj) < at{^kd{xk,S''x,) 



The fact that the support of {ai)i^z^ is contained in the balanced, convex and 
bounded set fl together with the recursion a["'' = J^j o,i-2jaf^ ^'^ (which amounts 
to the Chapman-Kolmogorov equations ([5])) implies supp(a(")) C (2" — l)il, see 
also 0]. 

Since the linear subdivision scheme with mask (a.i).igzs converges, we find a 
refinable function (/s : ^ R satisfying (fTT|) and , cf . Theorem [S] Recall that 
one obtains this refinable function as the limit of the linear scheme acting on the 
input data yj — Sjo , cf. Theorem [S] : 

sup|af)-^(z/2")| =e„ 0. (17) 
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Accordingly, setting A = {(fc,£) e x | max(p(i - 2"fc),/9(j - 2"^)) < 2" - 1}, 
we obtain 

{k,l)eA 



Ak,t)eA 



(18) 



\{k.e)eA 

Now, if j, fc, ^ e Z" are such that p(i - j) < 2, p(i - 2"fc) < 2" - 1 and p{j - 2"£) < 
2" — 1, one concludes 

p(fc -i)< ^{pii - 2"fc) + p(z - J-) + p(j - 

2; (19) 
< -(2(2" -1) + 2) = 2. 

Define 

■0(s, i) = X! "^(^ ~ ~ 

Then, since the refinable function is uniformly continuous, the property (|12p implies 
that for n large enough, 

a„ = inf V(s,i)>e>0- (20) 

p(t-s)<2-" + i 

By boundedness of H. we also obtain 

A/ = sup |z"' n (t + ri)| < 00. (21) 

Combining (HI]) with (HH) through (HH) further gives 

L»(S'"a;)= sup d(S'"a;,, S-^Xj) < 7„i:i(a;), (22) 

p(j-i)<2 

where 7„ = (1 — + 2e„ + APe^). Clearly, for uq large enough, 7 = 7„q < 1. 
Moreover, the estimate (|22p is uniform in a; (and even d). The same argument 
leading to the first inequality in (jl8D together with (1191) provides 

i:'(5"a;) < i:>(S''=a;) for m > k. 

Thus, for n e N one concludes: 

D{S"x) < jD{S'''""»x) 

< ^["/"o]£)(5"~"o["/"o]^ 

< 7[«/"olD(a;), 

which completes the proof. □ 

Recall that the tensor product (a (g) 6)jgxs+tof two masks {ai)i^z,= and {hj)j^jt is 
defined by 

(a ® 6)(ij-) = a,; • &j. 
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Lemma 10 ([5j). Suppose S and the corresponding D are as in Lemma\^ Define 
{bi)i^z "via bo = 1, 6_i = bi = ^, and 6^ = for \i\ > 1. Let T denote the 
barycentric scheme associated to the s-fold tensor power 6 (8) • ■ • 6. Then T is 
Lipschitz in the sense of ([5]) and converges on any Hadamard space. Moreover, 
there is C > Q such that doo{Sx,Tx) < C ■ D{x) for all x e ^°°(Z^X). 

Proof of Theorem [IJ Suppose S denotes the barycentric scheme associated to the 
nonnegative mask (ai)igz=- Under the assumption that the hnear counterpart of S 
converges, combining Lemmas 1^ and [TUl we find a function D : £°°{'E'^,X) — > R>o, 
a convergent scheme T : £°°(Z^X) £°°{Z'',X), and constants 7 < f and C > 
such that 

(i) There is a positive integer no such that D{S''^x) < 7["/"°lZ3(a;) for x 6 
^°°(Z^X) and n e N. 

(ii) T e Lipi(£°°(Z",X)) is convergent. 

(iii) d^{Sx,Tx) < C ■ D{x) for x G £°°(Z^ X). 

Thus, by Theorem [6j the scheme S converges. □ 

A well-known result from the linear theory is the following: 

Proposition 11. The univariate and linear scheme S associated to the mask 
(ai)iGZ converges if and only if there is 7 < 1 and C > such that 

sup \S"xi - S"'x,+i I < C • 7" sup \xi - Xi+i I for all n G Nq. 

Theorems [T] and IH] along with Lemma [H] put us in a position to generalize this 
statement to the setting of Hadamard spaces. Still need an easy auxiliary result. 

Lemma 12. Suppose {X, d) is a metric space, and let 

Doo{x) = sup d{xi,Xj). 

iK-j||oc<l 

Then a refinement scheme S is weakly contractive with respect to an admissible 
contractivity function if and only if there is 7 < 1 and C > such that 

DooiS'^x) < Cj^'D^ix) 

Proof. Suppose S is weakly contractive with respect to Dq , meaning there is ng G N 
and 7 < 1 such that £» o 5" < 7("/"o)£), is not difficult to see (cf. [5j) that there 
are r,R > such that 

rDn < L'oo < RDn- 

Observe that, since [n/no] > n/no - 1 one has 7["/"ol < ^^"/"o-i _ C7", where 
C = 7-1 and 7 = 7I/"" < 1. Moreover define C = Then 

Z?oo o 5" < RDn o 5" < i?7["/"°lz?n 

< -Cj'^D^ = C7"I?oo. 
r 

Now assume there is 7 < 1 and C > such that 

D^oS'' < Cj'^D^. 

Choose no G N such that < 1- From 1^ it follows that Dn o 5" < Dn- On 

the other hand, for n > uq we have [u/uq] < n — hq and thus 

Dn o 5" < -D^ o ^" < -Y'D^ 
r r 



□ 
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We are now able to generalize Proposition [TT] 

Theorem 7. The refinement scheme associated to {ai)i^z^ converges on arbitrary 
Hadamard spaces if and only if there is C > and 7 < 1 such that for all (X, d) 
Hadamard 

D^{S'"-x) <C ■-f''Doo{x) for allxee°°{I.',X), 
where, as above, Doo{x) = sup||,;_j||^<]^ (i(a;i, Xj). 

Proof. This follows from combining Lemma [T^ with Lemma [S] and Theorem [Sj □ 
We conclude this section with an approximation result for Lipschitz functions: 

Theorem 8. Suppose f : (R'*, || • ||) {X,d) is Lipschitz- continuous with constant 
C > 0, and S is a convergent barycentric scheme whose mask is supported on 
{a; G M** I ||x|| < r}. Sample f on the grid /iZ^, h > 0, via Xi = f{hi). Then 

doo{S^x{h-^-),f{-))<rC-h. 

Proof. Suppose n G Nq and i ^11 . Then by Lemma |S] and Theorem 2] we obtain 

d(5"x2.-.„/(/jz/2'=)) < a[lU^_^„^d{x,J{M/2'^)) 

||2n-fc,;_2'ij||<(2"-l)r 

< sup d{f{hj)J{hi/2'^)) 

\\i/2''-j\\<{l-2-")r 

<rC-h, 

from which the claim follows. □ 



3. LP-CONVERGENCE OF THE CHARACTERISTIC MARKOV CHAIN 

This short section clarifies the relationship between the stochastic convergence 
of the Markov chain associated to (ai)igzs and its counterpart in the theory of 
barycentric subdivision schemes. 

Lemma 13. Suppose supp(a) C CC\H , where C is a convex, balanced, and compact 
set. Let p : ^ M>o denote the Minkowski functional of C . Recall the notation 
Pi for the probability measure on (Z*)'*'° induced by the transition kernel P'^ and 
the initial distribution S^iy . Then 

p(i) < 2" =^ P,(x: e 2C) - 1. 

In other words, the Markov chain with deterministic initial condition Xq — i reaches 
2C within [log2(/o(i))] + 1 steps and remains in this set thereafter. 

Proof. Recall that since supp(a) C C, for any j e we obtain 

«!"Uy^O ^ p(*-2"j)<2". 
Thus the fact that p(i)/2" < 1 renders the right hand side of 



P.(X;: e \ 2C) = at\ 



pO)>2 



an empty sum, since p{i — 2'^j) < 2" implies 

p(j)<p(*-2'V)/2" + p(z)/2"<2. 



□ 
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Theorem 9. Let p e [l,oo). Suppose the characteristic Markov chain of 
(ai)igzs with deterministic initial condition £ g Z'' possesses a stationary distribu- 
tion TT in the sense that for all j g Z^, Wi^2"j ^ ^0 as n ~> oo. Then 
converges in LP(r2, P^; M*) if and only if there is k £ 1,'^ such that tt — Sk- In this 
case, 

Ei{\\X^ - k\\P) ^ asn^oo. 

Proof. Let p denote the Minkowski functional of a balanced, closed and convex set 
containing supp(a). Moreover, for £ N define 

An = eZ'xZ'l max(p(£ - 2^j),p{j - 2^i)) < 2" - 1}. 

Note that {i,j) G A„ implies that p{j) < 1 + '^^-^'ir"'" as well as p{i) < 1 + ^^^a^"*" . 
Thus there is a bounded set B such that 

U An C B. 

Moreover we have 

= ||^~JrP,(Xf„=^Kt = J)PKX,t=J) 



(23) 



Certainly, since B is bounded, the sequence 



e„ := sup (Itt, - a'f\,^\ + \tTj - a^ls-jl) 

converges to zero as n ^ oo. Consequently, we obtain Eg{\\X2n — > Cn, 

where 



E IK-Jr(«S". + «tU) (24) 

Thus, whenever there are integers i j such that tt^ > and tt^- > 0, Equation 
implies that iJ^dlXf^ — X^\\) is bounded away from zero asymptotically. Hence for 
LP-convergence of X^ we need the existence of some A: G with tt^ = Ski- 

Conversely, assume that tt = S^k}- Then since X^ kin distribution, X^ — > fc in 
probability. From Lemma [T51 we conclude that there is M > such that \\X^\\ < M 
holds Pf-almost surely for n E N. Hence for (5 > 0, 

ESK - k\n < I \\K - kVdPi + I \\K ~ fcfrfP^ 

J{||X^-fe||>(5} H\\X^-k\\<S} 

< (M + fc)PP,(||X--fc|| ><5)+<5P, 
showing that S^dlX,"^ ^ ^11) converges to zero. □ 
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Suppose now that the subdivision scheme associated to {ai)i^x' converges. Then 
a refinable function satisfying (fTTjl and exists. Substituting i for t in 

j 

and exploiting the fact that ip(i) = 1, we observe that tt^ = is a stationary 

distribution for X^. Moreover recall that a convergent scheme is called interpolatory 
if and only if there is fc e Z'' such that for j E 1/ , (p(j) = 5kj- Now Theorem [3] 
translates to the language of refinement schemes as follows: 

Corollary 14. Suppose the linear scheme associated to {ai)i^z^ converges, and 
p G [l,oo). Then the characteristic Markov chain of with deterministic initial 
condition ^ g converges in L^'(r2, P^; M*) if and only if the scheme is interpola- 
tory. In this case the L^-limit is a constant lattice point. 
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